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To provide a reliable and comprehensive data reference for core geometry design of graphite-moderated 
and low-enriched uranium fueled molten salt reactor, the influences of geometric parameters on temperature 
coefficient of reactivity (TCR) at an assembly level are characterized. The four-factor formula is introduced to 
explain how different reactivity coefficients behave in terms of fuel salt volume fraction and assembly size. The 
results show that fuel salt temperature coefficient (FSTC) is always negative due to a more negative fuel salt 
density coefficient in the over-moderated region or a more negative Doppler coefficient in the under-moderated 
region. Depending on the fuel salt channel spacing, the graphite moderator temperature coefficient (MTC) can 
be negative or positive. Further, an assembly with a smaller fuel salt channel spacing are more likely to exhibit 
a negative MTC. As fuel salt volume fraction increases, the negative FSTC weakens first and then increases, 
owing to the fuel salt density effect gradually weakening from negative feedback to positive feedback and then 
decreasing. Meanwhile, MTC weakens as the thermal utilization coefficient caused by the graphite temperature 
effect deteriorates. Thus, the negative TCR weakens first and then strengthens mainly because of the change 
in fuel salt density coefficient. As assembly size increases, the magnitude of FSTC decreases monotonously 
due to a monotonously weakened fuel salt Doppler coefficient, whereas MTC changes from gradual weakened 
negative feedback to gradual enhanced positive feedback. And then, the negative TCR weakens. Therefore, 
to achieve a proper negative TCR, particularly a negative MTC, an assembly with a smaller fuel salt channel 


spacing in the under-moderated region is strongly recommended. 
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I. INTRODUCTION 


Molten salt reactor (MSR) is the only liquid-fueled reactor 
in the generation IV power systems [1]. Its core structure and 
operation mode are obviously distinct from traditional pres- 
surized water reactor (PWR). A typical PWR uses water as 
coolant and moderator. As the core temperature of a PWR 
rises, the water density decreases and the ratio of water to 
uranium declines, but the fuel density remains almost con- 
stant. And PWR is usually designed in the under-moderated 
region to ensure a sufficiently negative moderator temperature 
coefficient [2]. Whereas a typical graphite-moderated MSR 
uses flowing fuel salt as coolant and solid graphite as mod- 
erator. As core temperature of an MSR rises, the density of 
graphite moderator remains essentially unchanged while the 
density of liquid fuel salt reduces, resulting in a higher ratio 
of graphite to nuclear fuel [3] and then significantly effects on 
temperature coefficient of reactivity (TCR). To make a reactor 
self-stable, it is crucial to maintain a proper negative temper- 
ature coefficient of reactivity (TCR) by adopting appropriate 
design parameters. 

The TCR of a graphite-moderated and liquid-fueled 
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MSR is typically divided into fuel salt temperature coeffi- 
cient (FSTC) and graphite moderator temperature coefficient 
(MTC). FSTC is usually negative, which can be further di- 
vided into fuel salt Doppler coefficient and fuel salt den- 
sity coefficient. However, MTC can be negative or posi- 
tive [4], and the latter poses a potential safety risk to reactor 
operation. Currently, most of liquid-fueled MSRs, such as 
MSRE [5], MSBR [6], and FUJI [7], were all designed in the 
under-moderated region. Nevertheless, these designs didn’t 
provided a systematic analysis to explain why the under- 
moderated region was chosen. 

Some research activities have been conducted on the in- 
fluencing factors of TCR in liquid-fueled MSRs, focusing on 
geometric parameters and fuel salt compositions. The geo- 
metric parameters [8—10] include fuel salt channel radius, fuel 
salt fraction, and lattice pitch, among others. The molar ratio 
of heavy metal (HM), the type of carrier salt, and the uranium 
enrichment are all important parameters associated with fuel 
salt compositions. The effects of fuel salt compositions on 
TCR for MSR have been studied using the six-factor formula 
[11-13]. Further research is however necessary to determine 
how geometric parameters affect TCR in an MSR core. 

Based on a graphite-moderated and low-enriched uranium 
fueled MSR, this work aims to provide a more responsible 
contribution to variations of various reactivity coefficients 
with geometric parameters including fuel salt volume frac- 
tion and assembly size. This paper is organized as follows. 
Section II provides a brief overview of calculation model and 
analysis methodology. Section III discusses the behaviors 
of various reactivity coefficients from the perspective of the 
four-factor formula. Section IV contains the concluded re- 
marks on TCR for liquid fueled MSRs. 


54 Il. CALCULATION MODEL AND ANALYSIS 
55 METHODOLOGY 
56 A. Calculation model 


s7 A graphite-moderated fuel salt assembly includes two re- 
ss gions (Figure 1): a central circular region filled with fuel salt 
sə (defined as fuel salt channel) and a hexagonal outer region 
eo composed of graphite moderator. The assembly geometry can 
61 be characterized by the side length of assembly (L), the op- 
62 posite side distance (P), the radius of fuel salt channel (R), 
63 and the fuel salt channel spacing (D) representing the dis- 
e4 tance between two adjacent fuel salt channels’ outer margins. 
65 It should be noted that the maximum fuel salt volume fraction 
66 (VF) is set to 50%. This is primarily because there is almost 
67 no Maxwell spectrum and MTC approaches zero when the 
e fuel salt volume fraction exceeds 50%. 

ə The adopted fuel salt composition in this work is LiF- 
70 BeF2-UF, (68-20-12 mole%). Previous studies have revealed 
71 more detailed analyses on the effect of the concentration of 
72 heavy metal on TCR [11, 12], which will not be discussed in 
73 this work. To reduce harmful neutron absorption of Li-6, the 
74 abundance of Li-7 is set to 99.995%. Low enriched uranium 
75 (<20% U-235 enrichment)) is used in this study for nonpro- 
ze liferation and availability. 


Fig. 1. Schematic of a graphite-moderated fuel salt assembly. 


B. Analysis methodology 


77 


7s This study was conducted for a single graphite-moderated 
73 fuel salt assembly that was infinitely reflected in both the 
so radial and axial directions. In this case, the moderator is 
s1 graphite, and the nuclear fuel is melted in the carrier salt [14], 
s2 Which includes F, Li, and Be, exhibits certain moderating 
83 properties. The Four-factor formula for an infinite medium 
84 [15], kine=epnf, is used to qualitatively understand the mecha- 
es nisms of TCR of an MSR at assembly level. Here, £, p, 7 and 
ss f represent the fast fission factor, the resonance escape prob- 
s7 ability, the thermal reproduction factor, and the thermal uti- 
ss lization factor, respectively. The four factors and their corre- 
s sponding reactivity coefficients for a graphite-moderated fuel 
90 Salt assembly with low enriched uranium are detailed as fol- 
ə lows. It is crucial to keep in mind that the nuclear fuels (or 
ve heavy metals) are integrated into the fuel salt, the fission or 


ə absorption contributions of them should be subtracted from 
o4 the fuel salt when calculating the four factors, especially for 
ə the fast fission factor and the thermal reproduction factor. 


æ The fast fission factor, € = FM/F™, is defined as the 
97 ratio of the number of neutrons produced by fissions at all 
ə energies (total production) to that produced by thermal fis- 
ə Sion (thermal production). The thermal reproduction fac- 
tor, n = FM/AM, is the ratio of the number of neutrons 
produced by thermal fission (thermal production) to that of 
thermal neutrons absorbed in nuclear fuel (thermal absorp- 
tion). Here, F}, Ff, and A, represent the total production, the 
thermal production and the thermal absorption, respectively. 
Since fast fission occurs primarily in U-238 but also in U- 
235 when low enriched uranium is used, the subscript "HM" 
includes U-235 and U-238. 

The resonance escape probability, p = Att / Att = Ast + 
A‘, is described as the ratio of neutrons reaching thermal 
energies to fast neutrons slowing down. The thermal utiliza- 
tion factor, f = AIM/ Att = AM /( AEM. ĄFliBe 1 Af?) — 
AHM /(AHM 4 ATMBes) is represented as the ratio of ther- 
mal neutrons absorbed in heavy metal to those absorbed in 
all the materials in the assembly. Here, A; represents the to- 
tal absorption and A‘ equals 1.0. The superscripts "HM", 
"FLiBe", "gra" and "tot" represent the thermal absorption of 
heavy metal, carrier salt, graphite and all the materials in the 
assembly, respectively. 


The temperature coefficient of reactivity, œr, defined as the 
change in reactivity caused by the change in temperature of 
all the materials in the assembly [16], is expressed by Eq. (1) 
when king approaches 1.0. Temperature (T) rises in 100 K 
increments from 780 K to 1080 K considering the melting 
124 point of the fuel salt and the possible temperature range of 
125 MSR. 


1 dking 
kint dT 


126 QT (1) 


127 According to Eq. (2), each TCR can be decomposed into four 
128 reactivity coefficients [17], including fast fission coefficient 
129 (af), resonance escape coefficient (at), thermal reproduction 


130 Coefficient (at), and thermal utilization coefficient (ad), 


(2) 
132 Eqs. (3)~ (6) presents the calculation methods for four reac- 
133 tivity coefficients. Here, "d" denotes the absolute change in 
134 parameter caused by temperature change and "A" denotes the 
135 rate at which the parameter changes with temperature. 


a1 ar = ah + ak +a} + af 


1 de HM HM 1 de 1 dk 

136 a= edt = AFP’-AFy" = FEM F FM dT 6) 
1dP_ aas 1 daf" dA 

r= ——_= AAS AAP = ag? ap 

137 QT pdt 1 + 1 Asat qar. dT dT ) ( ) 
1d HM HM 1 aM 1 dai 
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It’s important to note that U-235 plays a dominate role in 
feedback, particularly for af and ał, because its microscopic 
cross section is significantly higher than that of U-238 over 
the thermal energy region, particularly when neutron energy 
is less than 0.1 eV. Therefore, in the following discussions, 
we will focus on U-235’s contribution to the fast fission coef- 
ficient (af) and the thermal reproduction coefficient (a). 

MCNP5S handles neutronic calculations such as critical- 
ity and reactivity coefficients. The FMn tally multiplier 
card is used to calculate cross sections and factors in the 
four-factor formula. To perform accurate calculations with 
uranium-based fuel, a compact ENDF (ACE) format cross 
section library with continuous energy was chosen based on 
the ENDF/B-VII library. In order to increase computational 
accuracy and efficiency, each criticality calculation is ex- 
cused to skip 50 cycles and run a total of 1050 cycles with 
100,000 neutrons per cycle. The maximum computing time 
for one criticality calculation using 12 processors is less than 
36 hours. The estimated standard deviation of kint is about 5 
pem. 


HI. RESULTS AND DISCUSSIONS 


A. Critical parameters 


Considering the characteristics of online refueling and re- 
processing for a liquid-fueled MSR, the initial excess reactiv- 
ity of a graphite-moderated fuel salt assembly can be set low 
(King =1.02~1.03). And then, the required critical enrichment 
of U-235 is searched by varying assembly size (L) from 3.0 
cm to 24.0 cm and fuel salt volume fraction (VF) from 1% to 
50%. The fuel salt channel spacing is closely related to the 
required critical enrichment of U-235 and the temperature co- 
efficient of reactivity, and its variation with fuel salt volume 
fraction and assembly size is shown in Figure 2 (a). Increas- 
ing fuel salt volume fraction and decreasing assembly size 
can result in a decrease in fuel salt channel spacing. Mean- 
while, the larger the assembly size, the faster fuel salt channel 
spacing decreases with increasing fuel salt volume fraction. 

The variations of the infinite multiplication factor at a fixed 
U-235 enrichment (19.75 wt%) and the required critical en- 
richment of U-235 with fuel salt volume fraction and assem- 
bly size are given in Figure 2 (b) and Figure 2 (c), respec- 
tively. When fuel salt volume fraction is raised for a fixed 
assembly size, the neutron spectrum hardens, and the resi- 
dence time of neutrons in high energy region of fast fission 
is prolonged. Hence, the fast neutron multiplication effect 
is enhanced, and the fast fission factor increases. Mean- 
while, as graphite moderator volume fraction decreases, the 
ratio of graphite to heavy metal declines, and the probability 
of neutrons being absorbed by nuclear fuel or moderator in 
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Fig. 2. Variations of fuel salt channel spacing (a), kint (b), critical 
enrichment of U-235 (c) and EALF (d) with fuel salt volume fraction 
for different assembly sizes. 


epi-thermal neutron energy region increases, making it more 
difficult for neutrons to escape from resonance region, and 
thus the resonance escape probability decreases. The thermal 
reproduction factor, on the other hand, varies only slightly 
with fuel salt volume fraction, with a maximum variation of 
less than 0.6%. As the ratio of fuel salt volume to graphite 
volume further increases, the thermal absorption of fuel salt 
exceeds that of graphite, and the thermal utilization factor 
rises. Based on the variations of four factors, the infinite mul- 
tiplication factor first increases corresponding to the under- 
moderated region and then decreases corresponding to the 
over-moderated region as fuel salt volume fraction increases. 
On the contrary, the required critical enrichment of U-235 de- 
creases in the over-moderated region and then increases in the 
under-moderated region. 


For a constant fuel salt volume fraction, as assembly size 
increases, the graphite through which the incident neutrons 
thickens, neutrons can be slowed down even better, neutron 
spectrum softens. In this case, the fast neutron multiplication 
effect weakens, and the fast fission factor decreases. Mean- 
while, the absorption probability of resonance neutrons gen- 
erated in graphite by heavy metal decreases, resulting in an 
increase in resonance escape probability [18]. The thermal 
reproduction factor is also not visibly affected by the assem- 
bly size, with a maximum variation of less than 0.4%. As 
fuel salt channel spacing increases with increasing assembly 
size, the thermal neutron absorption of nuclear fuel decreases 
while that of graphite moderator increases, and then the ther- 
mal utilization factor decreases. As assembly size increases, 
a slight decrement in the infinite multiplication factor and a 
slight increment in the required critical enrichment of U-235 
are displayed in the over-moderated region. Meanwhile, the 
infinite multiplication factor rises while the required critical 
enrichment of U-235 falls in the under-moderated region. 


Neutron spectrum is essential for determining the influ- 
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ence factors of TCR. To quantitatively characterize the neu- 
tron spectrum in a reactor, a spectrum factor defined as the 
energy corresponding to the average neutron lethargy caus- 
ing fission (EALF) [12] is introduced and its variation with 
fuel salt volume fraction and assembly size is displayed in 
Figure 2 (d). As fuel salt volume fraction increases for a con- 
stant assembly size, the required critical enrichment of U-235 
increases especially in the under-moderated region, allowing 
for more production of fast neutrons. As fuel salt channel 
spacing decreases, the fast neutrons released from fuel salt 
cannot be fully moderated. These two factors harden neu- 
tron spectrum, leading to an increase in EALF. The effect of 
assembly size on EALF is associated with the moderated re- 
gion. First, when the assembly tends to the over-moderated 
region, the graphite parasitic absorption becomes stronger, 
the likelihood of fast neutrons being slowed down to ther- 
mal neutrons decreases. And, a slightly increasing required 
critical enrichment of U-235 results in more fast neutron gen- 
eration as assembly size increases. Thus, neutron spectrum 
hardens and EALF increases as assembly size increases in the 
over-moderated region. Second, when the assembly tends to 
the under-moderated region, the graphite dominates the scat- 
tering reaction. As assembly size increases, the likelihood 
of fast neutrons colliding with graphite nuclide increases, as 
does the possibility of fast neutrons being slowed down into 
thermal neutrons. At the same time, an increasing assembly 
size causes a decrease in the required critical enrichment of 
U-235 and less fast neutron generation. Thus, neutron spec- 
trum softens and EALF decreases as assembly size increases 
in the under-moderated region. 


B. Fuel salt Doppler coefficient 


The Doppler effect, which can be used to explain the 
Doppler broadening of the resonance capture cross-sections 
of nuclear fuel since it determines the fuel salt temperature 
coefficient. From the point of view of the four-factor formula, 
the fuel salt Doppler coefficient can be decomposed into four 
reactivity coefficients, i.e., fast fission coefficient, resonance 
escape coefficient, thermal reproduction coefficient, and ther- 
mal utilization coefficient, as demonstrated in Eq. (2). Their 
variations with fuel salt volume fraction and assembly size 
are presented in Figure 3. 

The Doppler broadening of resonance caused by an in- 
crease in fuel salt temperature leads to a hardening neutron 
spectrum, which causes the residence time of neutrons in 
high energy region of fast fission to be longer, the fast neu- 
tron multiplication effect to be enhanced, and the fast fission 
factor increases. Thus, the fast fission coefficient in Doppler 
coefficient, a7(dop), is positive. Its magnitude is primarily 
determined by the difference between U-235’s thermal pro- 
duction and U-235’s total production, and it is very closely 
related to the change of neutron spectrum. When fuel salt 
volume fraction or assembly size changes, if neutron spec- 
trum becomes harder and EALF increases, which means that 
the proportion of U-235’s thermal production in U-235’s total 
production decreases, and then the gap between the variations 
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Fig. 3. Variations of different reactivity coefficients caused by the 
Doppler effect with fuel salt volume fraction and assembly size. 
(a) at(dop): fast fission coefficient in Doppler coefficient; (b) 
aĵ (dop): resonance escape coefficient in Doppler coefficient; (c) 
ai(dop): thermal reproduction coefficient in Doppler coefficient; 
(d) af (dop): thermal utilization coefficient in Doppler coefficient; 
and (e) ar(dop): Doppler coefficient. 


in U-235’s total production and U-235’s thermal production 
becomes wider. It implies that the magnitude of fast fission 
coefficient in Doppler coefficient grows. By contrary, when 
neutron spectrum softens, the magnitude of fast fission co- 
efficient in Doppler coefficient decreases. As a result, Fig- 
ure 3 (a) shows that for a constant assembly size, the magni- 
tude of fast fission coefficient in Doppler coefficient increases 
monotonously as fuel salt volume fraction increases. As as- 
sembly size increases for a fixed fuel salt volume fraction, the 
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magnitude of fast fission coefficient in Doppler coefficient in- 
creases slightly in the over-moderated region but significantly 
decreases in the under-moderated region. Overall, the magni- 
tude of fast fission coefficient in Doppler coefficient is posi- 
tively correlated with EALF. 


When fuel salt temperature rises, the Doppler broadening 
effect enhances the resonance absorption of nuclear fuel and 
lowers the likelihood of neutrons passing through the epi- 
thermal region to thermal region, the thermal absorption of 
fuel salt decreases. Nevertheless, the thermal absorption of 
graphite may decrease for a smaller fuel salt channel spacing 
or may increase for a larger fuel salt channel spacing. This 
is because the thermal neutrons generated mainly in graphite 
have higher probability to be absorbed by graphite itself when 
fuel salt channel spacing is relatively large. Since the Doppler 
effect mainly affects the heavy metals in fuel salt, the change 
in thermal absorption of fuel salt is usually at least one or- 
der of magnitude greater than the change in thermal absorp- 
tion of graphite. Thus, the sum of changes in thermal ab- 
sorptions of fuel salt and graphite is always negative. And 
consequently, the resonance escape coefficient in Doppler co- 
efficient, at (dop), is negative. From Figure 3 (b), it can be 
concluded that the variation of resonance escape coefficient 
in Doppler coefficient comes mainly from the sum of the ther- 
mal absorptions of fuel salt and graphite, and to a lesser de- 
gree, from their changes. That is, for a constant assembly 
size, as fuel salt volume fraction increases, a harder neutron 
spectrum makes the sum of the thermal absorptions of fuel 
salt and graphite lower. Hence, a larger variation in resonance 
escape probability is revealed, and a stronger resonance es- 
cape coefficient in Doppler coefficient is presented for a larger 
fuel salt volume fraction corresponding to a smaller fuel salt 
channel spacing. As assembly size increases for a fixed fuel 
salt volume fraction, the fuel salt channel spacing widens, the 
sum of the thermal absorptions of fuel salt and graphite grad- 
ual increases, weakening the resonance escape coefficient in 
Doppler coefficient. 


The Doppler broadening effect hardens neutron spectrum, 
increasing the resonance absorption of nuclear fuel and de- 
creasing the probability of thermal neutrons being absorbed 
by nuclear fuel. Both U-235’s thermal production and U- 
235’s total absorption decrease. Since the microscopic fission 
cross section of U-235 is greater than its microscopic cap- 
ture cross section, especially in thermal neutron region, and 
the average number of fission neutrons of U-235 is generally 
greater than 2.0, the absolute decrement of thermal produc- 
tion is always stronger than that of thermal absorption. For 
most combinations of fuel salt volume fraction and assembly 
size, despite the fact that U-235’s thermal production is larger 
than U-235’s thermal absorption, the absolute variation in U- 
235’s thermal production is still larger than that in U-235’s 
thermal absorption, and then the thermal reproduction factor 
decreases and further leads to a negative thermal reproduc- 
tion coefficient in Doppler coefficient (a? (dop)). However, 
for an assembly with a very small fuel salt channel spacing 
(e.g., L=3 cm and VF =40%), the graphite thickness that neu- 
trons passes is relatively small, neutrons are not sufficiently 
slowed down, the probability of thermal neutron production 
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is reduced, as does U-235’s absorption (both fission and cap- 
ture). In this case, the difference between U-235’s thermal 
production and its thermal absorption is greater than the dif- 
ference between the decrement of U-235’s thermal produc- 
tion and that of U-235’s thermal absorption. Thus, the varia- 
tion in U-235’s thermal production has a lower absolute value 
than the variation in U-235’s thermal absorption, increasing 
the thermal reproduction factor and causing a positive ther- 
mal reproduction coefficient in Doppler coefficient. It can be 
concluded from Figure 3 (c) that the magnitude of thermal 
reproduction coefficient in Doppler coefficient increases first 
and then decreases as fuel salt volume fraction increases but 
increases monotonously with increasing assembly size. First, 
as fuel salt volume fraction increases, a hardening neutron 
spectrum causes decrements in both U-235’s thermal pro- 
duction and its thermal absorption, followed by increments 
in both the variation in U-235’s thermal production and that 
in U-235’s thermal absorption. When fuel salt volume frac- 
tion begins to rise, the proportion of thermal neutrons is rela- 
tively high, the increase of the variation in U-235’s thermal 
production is greater than that in U-235’s thermal absorp- 
tion, causing the magnitude of thermal reproduction coeffi- 
cient in Doppler coefficient to increase slightly. However, 
as fuel salt volume fraction continues to increase, the share 
of thermal neutron absorption decreases while the share of 
resonance neutron absorption and fast neutron absorption in- 
creases. At this point, the increase of the variation in U-235’s 
thermal absorption gradually exceeds that in U-235’s thermal 
production, causing the magnitude of thermal reproduction 
coefficient in Doppler coefficient to decrease. And the neg- 
ative thermal reproduction coefficient in Doppler coefficient 
strengths first and then weakens as fuel salt volume fraction 
increases for a constant assembly size. Second, the change in 
thermal reproduction coefficient in Doppler coefficient with 
assembly size is divided into two phases according to fuel 
salt volume fraction. When fuel salt volume fraction is very 
small (e.g., WF=1%) or very large (e.g., VF=40%), the ab- 
solute changes of U-235’s thermal production and U-235’s 
thermal absorption increase with increasing assembly size, 
but the former increases faster than the latter. As a result, 
the increment of the variation in U-235’s thermal production 
is greater than that in U-235’s thermal absorption, the mag- 
nitude of thermal reproduction coefficient in Doppler coef- 
ficient increases as assembly size increases. When fuel salt 
volume fraction is close to the optimal moderated zone (e.g., 
VF=15%), an increasing resonance absorption share causes 
the absolute changes of U-235’s thermal production and that 
of U-235’s thermal absorption to decrease with increasing 
assembly size, with the latter decreasing faster than the for- 
mer. And then, the decrease of the variation in U-235’s ther- 
mal absorption is greater than that in U-235’s thermal pro- 
duction, the magnitude of thermal reproduction coefficient 
in Doppler coefficient is still increasing as assembly size in- 
creases. Therefore, when fuel salt volume fraction is constant, 
the magnitude of thermal reproduction coefficient in Doppler 
coefficient increases monotonously with increasing assembly 
size. It should be noted that the magnitudes of thermal repro- 
duction coefficient in Doppler coefficient for all combinations 
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of fuel salt volume fraction and assembly size are quite small 
(< 0.7 pem/K). 


The thermal absorption of fuel salt decreases as fuel salt 
temperature rises owing to a hardening neutron spectrum 
caused by the Doppler effect. After that, the thermal absorp- 
tions of heavy metal and carrier salt decrease. The variation 
in thermal absorption of graphite is related to the fuel salt 
channel spacing. First, for assembly with a smaller fuel salt 
channel spacing, the thermal absorption of graphite is primar- 
ily affected by the hardening neutron spectrum, and the ther- 
mal absorption of graphite decreases, as does the sum of ther- 
mal absorptions of carrier salt and graphite. Second, as fuel 
salt channel spacing increases, neutrons have more chances 
of colliding with graphite nuclei, the probability of neutrons 
slowing down into thermal neutrons increases, and the ther- 
mal absorption of graphite increases. In this case, the sum of 
the thermal absorptions of carrier salt and graphite may in- 
crease in the over-moderated region or decrease in the under- 
moderated region. In either case, however, the change in ther- 
mal absorption of heavy metal is always stronger than that in 
the sum of thermal absorptions of carrier salt and graphite 
caused by the Doppler effect, resulting in a decrease in ther- 
mal utilization factor and a negative thermal utilization coef- 
ficient in Doppler coefficient (at (dop)). The magnitude of 
thermal utilization coefficient in Doppler coefficient is pri- 
marily related to the ratio of the sum of the thermal absorp- 
tions of carrier salt and graphite to the thermal absorptions 
of all materials according to Eq. 6. From Figure 3 (d), for 
a fixed assembly size, as fuel salt volume fraction increases, 
fuel salt channel spacing decreases, the graphite volume frac- 
tion decreases, and the ratio of the sum of the thermal absorp- 
tions of carrier salt and graphite to the thermal absorptions 
of all materials is less than 1.0 and becomes smaller, while 
the negative thermal utilization coefficient in Doppler coef- 
ficient weakens. Meanwhile, for a constant fuel salt volume 
fraction, as assembly size increases, the graphite thickness 
through which neutrons passed between two collisions with 
nuclear fuel increases, implying that more thermal neutrons 
may be absorbed by graphite. Then, the ratio of the sum of 
the thermal absorptions of carrier salt and graphite to the ther- 
mal absorptions of all materials grows larger and larger, and 
the negative thermal utilization coefficient in Doppler coeffi- 
cient strengths. To summarize, the negative thermal utiliza- 
tion coefficient in Doppler coefficient is proportional to fuel 
salt channel spacing. 


Figure 3 (e) depicts the variation of Doppler coefficient 
with fuel salt volume fraction and assembly size. The be- 
havior of Doppler coefficient is primarily governed by the 
variations of fast fission coefficient, resonance escape coef- 
ficient, and thermal utilization coefficient. Because the sum 
of the negative resonance escape coefficient and the negative 
thermal utilization coefficient is always greater than the pos- 
itive fast fission coefficient, the Doppler coefficient is always 
negative. The magnitude of Doppler coefficient varies with 
fuel salt volume fraction, mainly depending on the compe- 
tition between fast fission coefficient and resonance escape 
coefficient. Increasing fuel salt volume fraction increases the 
magnitude of Doppler coefficient firstly due mainly to an en- 
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hancement in resonance escape coefficient, and it reaches a 
maximum negative value before beginning to decrease with 
further increases in fuel salt volume fraction due to an en- 
hancing positive fast fission coefficient. As assembly size 
increases, the magnitude of fuel salt Doppler coefficient in- 
creases in the over-moderated region due to a slight enhance- 
ment in thermal utilization coefficient and decreases in the 
under-moderated region due to a gradually weakened reso- 
nance escape coefficient. 


C. Fuel salt density coefficient 


The fuel salt density effect occurs because a small amount 
of fuel salt is ejected from the reactor core as its density de- 
creases with increasing fuel salt temperature. The decrease 
in fuel salt density results in two major effects: the first is 
an increase in collision probability between fast neutron and 
graphite nuclei as a result of the reduced collision probabil- 
ity between fast neutron and nuclear fuel, and the second is a 
reduction in resonance absorption of nuclear fuel, which re- 
sults in more fast neutrons being slowed down to thermal neu- 
trons. Both effects can soften neutron spectrum and influence 
fuel salt density coefficient. Similar to fuel salt Doppler co- 
efficient, fuel salt density coefficient can also be divided into 
four parts, namely fast fission coefficient, resonance escape 
coefficient, thermal reproduction coefficient, and thermal uti- 
lization coefficient caused by the fuel salt density effect. The 
variations of fuel salt density coefficient and its four reactivity 
coefficients are presented in Figure 4. 

When fuel salt temperature increases, a softening neutron 
spectrum caused by the fuel salt density effect weakens the 
fast neutron multiplication effect, and then the fast fission fac- 
tor reduces. Thus, the fast fission coefficient in density coeffi- 
cient, a7(den), is negative. The variation of fast fission coeffi- 
cient in density coefficient with fuel salt volume fraction and 
assembly size is similar to that in Doppler coefficient, and it 
is primarily determined by the shift of neutron spectrum. As 
shown in Figure 4 (a), for a constant assembly size, as fuel 
salt volume fraction increases, neutron spectrum hardens, the 
proportion of thermal production in total production for U- 
235 decreases, the gap between the variation in U-235’s total 
production and that in U-235’s thermal production enlarges, 
and then the fast fission coefficient in density coefficient is en- 
hanced. Similarly, as assembly size increases, the magnitude 
of fast fission coefficient in density coefficient either increases 
because of a gradual hardening neutron spectrum in the over- 
moderated region or decreases because of a gradual softening 
neutron spectrum in the under-moderated region. In general, 
the magnitude of fast fission coefficient in density coefficient 
is positively correlated with EALF. 

As fuel salt temperature rises, a softening neutron spectrum 
caused by the fuel salt density effect increases the thermal ab- 
sorption of graphite. The change in thermal absorption of fuel 
salt is influenced by both density reduction and spectrum soft- 
ening, which is closely related to the moderated region. In the 
over-moderated region, the change in thermal absorption of 
fuel salt is due mainly to a decrease in fuel salt density, result- 
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Fig. 4. Variations of different reactivity coefficients caused by the 
fuel salt density effect with fuel salt volume fraction and assem- 
bly size. (a) at(den): fast fission coefficient in density coefficient; 
(b) af(den): resonance escape coefficient in density coefficient; (c) 
a}(den): thermal reproduction coefficient in density coefficient; (d) 
af (den): thermal utilization coefficient in density coefficient; and 
(e) ar(den): density coefficient. 


ing in a numerical decrease. In this case, a lower fuel salt vol- 
ume fraction corresponds to a higher graphite volume share, 
and the decrease of thermal absorption of fuel salt is smaller 
than the increase of thermal absorption of graphite. Thus, the 
sum of thermal absorptions of fuel salt and graphite is pos- 
itive. In the under-moderated region, the change of thermal 
absorption of fuel salt is primarily due to a softening neutron 
spectrum, and its value will increase, as does the sum of the 
thermal absorptions of fuel salt and graphite. The resonance 
escape probability increases in either case, and the resonance 
escape coefficient in density coefficient, af(den), is positive. 
From Figure 4 (b), the magnitude of resonance escape coef- 
ficient in density coefficient is also associated with fuel salt 
channel spacing, and its variation with fuel salt volume frac- 
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tion and assembly size is also determined by the sum of the 
thermal absorptions of fuel salt and graphite. As fuel salt vol- 
ume fraction increases or assembly size decreases, the sum 
of thermal absorptions of fuel salt and graphite increases, as 
does the magnitude of resonance escape coefficient in density 
coefficient. In short, the magnitude of resonance escape co- 
efficient in density coefficient is inversely related to the fuel 
salt channel spacing. 


The fuel salt density effect causes a decrease in fuel salt 
density and then a softening neutron spectrum as fuel salt 
temperature increases. These two variables can determine 
whether the thermal reproduction coefficient in density coef- 
ficient, a/(den), is positive or negative. First, the density re- 
duction is usually reflected in the region with a relatively soft 
neutron spectrum (such as VF=1%). In this case, even though 
the decrement of U-235’s thermal production is greater than 
that of U-235’s thermal absorption, the decrement of the vari- 
ation in U-235’s thermal production is smaller than that in 
U-235’s thermal absorption due to the fact that U-235’s ther- 
mal production is much greater than U-235’s thermal absorp- 
tion, and then a positive thermal reproduction coefficient in 
density coefficient is presented. Second, a softening neutron 
spectrum causes increases of both U-235’s thermal produc- 
tion and U-235’s thermal absorption, but the former is faster 
than the latter. The increase in the variation in U-235’s ther- 
mal production and that in U-235’s thermal absorption in this 
case is closely related to the fuel salt channel spacing. For 
a larger fuel salt channel spacing, the difference between U- 
235’s thermal production and U-235’s thermal absorption has 
no effect on the gap between the increment of U-235’s ther- 
mal production and that of U-235’s thermal production, and 
then the increment in the variation in U-235’s thermal pro- 
duction is still greater than that in U-235’s thermal absorp- 
tion, resulting in a positive thermal reproduction coefficient 
in density coefficient. For a smaller fuel salt channel spac- 
ing, the difference between U-235’s thermal production and 
U-235’s thermal absorption is greater than the difference be- 
tween their respective increments, and then the increment in 
the variation in U-235’s thermal production is smaller than 
that in U-235’s thermal absorption, revealing a negative ther- 
mal reproduction coefficient in density coefficient. With in- 
creasing fuel salt volume fraction, the positive thermal re- 
production coefficient caused by the density reduction effect 
gradually turns to the positive thermal reproduction coeffi- 
cient caused by the spectrum softening effect. Therefore, a 
negative feedback is attained. Figure 4 (c) shows that with 
increasing fuel salt volume fraction, the positive thermal re- 
production coefficient in density coefficient firstly strengths, 
then weakens and becomes negative. The variation of thermal 
reproduction coefficient in density coefficient with assembly 
size is related to the moderated region. There is no signif- 
icant difference in thermal reproduction coefficient in den- 
sity coefficient between different assembly sizes in the over- 
moderated region. An increasing assembly size in the under- 
moderated region results in a weakening negative feedback or 
a strengthening positive feedback. When assembly is near the 
optimal moderated region, the positive feedback caused by 
the softening spectrum effect gradually weakens as assembly 
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size increases. Overall, the thermal reproduction coefficient 
in density coefficient contributes very little (< 0.1 pem/K) to 
fuel salt density coefficient. 


As fuel salt temperature rises, its density decreases and 
neutron spectrum softens. Graphite’s thermal absorption in- 
creases due to a softening neutron spectrum, and then the ther- 
mal absorption of graphite grows larger. The moderated re- 
gions affect the change in thermal absorption of fuel salt. In 
the over-moderated region, the thermal absorption of fuel salt 
decreases mainly due to a reduction of fuel salt density, as will 
the thermal absorption of heavy metal and that of carrier salt. 
In this case, because the increase in thermal absorption of 
graphite is greater than the decrease in thermal absorption of 
carrier salt as a result of a larger graphite volume fraction, the 
sum of thermal absorptions of graphite and carrier salt shows 
an increase effect, and then the thermal utilization factor de- 
creases. In the under-moderated region, the thermal absorp- 
tion of heavy metal and that of carrier salt increase mainly 
due to a softening neutron spectrum, but the reduction in fuel 
salt density attenuates the increase rate of thermal absorption 
of heavy metal. Since the increase in thermal absorption of 
heavy metal is less than that in the sum of the thermal absorp- 
tions of carrier salt and graphite, the thermal utilization factor 
decreases. In any case, the thermal utilization coefficient in 
density coefficient, af (den), is always negative (Figure 4 (d)). 
Similarly, the magnitude of thermal utilization coefficient in 
density coefficient is closely related to the ratio of the sum of 
the thermal absorptions of carrier salt and graphite to the ther- 
mal absorptions of all materials. As fuel salt volume fraction 
increases for a constant assembly size, the ratio of the sum 
of the thermal absorptions of carrier salt and graphite to the 
thermal absorptions of all materials gradually decreases with 
a decreasing graphite volume fraction, and then the negative 
thermal utilization coefficient in density coefficient weakens. 
At a fixed fuel salt volume fraction, the difference in ther- 
mal utilization coefficient in density coefficient for different 
assembly sizes can be negligible. 


The change in fuel salt density coefficient with fuel salt 
volume fraction and assembly size is depicted in Figure 4 
(e). Because the magnitude of thermal reproduction coeffi- 
cient is very small, the feedback and the magnitude of density 
coefficient vary with fuel salt volume fraction and assembly 
size, owing primarily to variations in fast fission coefficient, 
resonance escape coefficient, and thermal utilization coeffi- 
cient. First, in the over-moderated region, because the feed- 
backs of fast fission coefficient and resonance escape coef- 
ficient are opposite and their magnitudes are close, the den- 
sity coefficient exhibits a negative feedback due to a more 
negative thermal utilization coefficient. Second, in the under- 
moderated region, the magnitude of thermal utilization coef- 
ficient is small in comparison to the magnitudes of fast fission 
coefficient and resonance escape coefficient, and the feed- 
back and the magnitude of density coefficient are primarily 
determined by the latter two. Because the negative fast fis- 
sion coefficient is weaker than the positive resonance escape 
coefficient in this case, the feedback of density coefficient is 
positive. Finally, the density coefficient is more sensitive to 
fuel salt volume fraction rather than assembly size. With in- 
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creasing fuel salt volume fraction, the negative density coeffi- 
cient decreases firstly due to a decreasing thermal utilization 
coefficient, then turns into an increasing positive density co- 
efficient due to an increasing positive resonance escape coef- 
ficient, and finally the positive density coefficient decreases 
due to an increasing negative fast fission coefficient. 


D. Fuel salt temperature coefficient 


Fuel salt temperature coefficient (FSTC) is a cumulative 
effect of fuel salt Doppler coefficient and fuel salt density co- 
efficient, which can also be divided into four parts, namely 
fast fission coefficient, resonance escape coefficient, thermal 
reproduction coefficient, and thermal utilization coefficient 
caused by the fuel salt temperature effect. Figure 5 shows the 
variations of four reactivity coefficients that are responsible 
for the behavior of FSTC. 

The change in fast fission factor caused by the temperature 
change of fuel salt, œẸ(salt), is a sum of fuel salt Doppler 
effect and fuel salt density effect. Since the contribution of 
the positive Doppler effect is smaller than the negative fuel 
salt density effect, the fast fission coefficient in FSTC is neg- 
ative. The variation of fast fission coefficient in FSTC with 
fuel salt volume fraction and assembly size is very similar to 
that in fuel salt density effect, and is displayed in Figure 5 (a). 
When neutron spectrum hardens (or softens) due to variations 
in fuel salt volume fraction or assembly size, the EALF raises 
(or falls), the magnitude of fast fission coefficient in FSTC 
increases (or decreases). That is, the magnitude of fast fission 
coefficient in FSTC correlates positively with EALF. 

The change in resonance escape probability caused by the 
temperature change in fuel salt, af (salt), is also a sum of fuel 
salt Doppler effect and fuel salt density effect. For most com- 
binations of fuel salt volume fraction and assembly size, the 
positive resonance escape coefficient caused by the fuel salt 
density effect is stronger than the negative resonance escape 
coefficient caused by the Doppler effect, resulting in a pos- 
itive resonance escape coefficient in FSTC. However, for an 
assembly with a smaller assembly size (e.g., L=3 cm) and a 
smaller fuel salt volume fraction (e.g., VF<10%), the posi- 
tive resonance escape coefficient caused by the fuel salt den- 
sity effect is weaker than the negative resonance escape coef- 
ficient caused by the Doppler effect, resulting in a relatively 
weaker negative resonance escape coefficient in FSTC. From 
Figure 5 (b), at a constant assembly size, the negative reso- 
nance escape coefficient in FSTC becomes weaker and turns 
to a more stronger positive feedback with increasing fuel salt 
volume fraction. This is mainly due to a gradual increased 
positive resonance escape coefficient caused by the fuel salt 
density effect. The magnitude of resonance escape coefficient 
in FSTC with a fixed fuel salt volume fraction is divided into 
two main regions by the change of assembly size. In the over- 
moderated region, as assembly size increases, the resonance 
escape coefficient in FSTC shifts from weaker negative feed- 
back to stronger positive feedback due to a weakening neg- 
ative resonance escape coefficient caused by the Doppler ef- 
fect. In the under-moderated region, a decreasing positive res- 
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Fig. 5. Variations of different reactivity coefficients caused by the 
fuel salt temperature effect as function of fuel salt volume fraction 
and assembly size. (a) œî(salt): fast fission coefficient in FSTC; (b) 
a#(salt): resonance escape coefficient in FSTC; (c) a#(salt): ther- 


mal reproduction coefficient in FSTC; (d) af (salt): thermal utiliza- 
tion coefficient in FSTC; and (e) œr(salt): FSTC. 


onance escape coefficient caused by the fuel salt density ef- 
fect causes the positive resonance escape coefficient in FSTC 
to gradually decrease with increasing assembly size. In gen- 
eral, the larger the EALF, the greater the possibility of a pos- 
itive resonance escape coefficient in FSTC and the stronger 
the positive resonance escape coefficient in FSTC. 


The thermal reproduction coefficient in FSTC is also a sum 
of fuel salt Doppler effect and fuel salt density effect. Because 
the magnitude of thermal reproduction coefficient caused by 
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the fuel salt density effect is very small, the feedback and 
the magnitude of thermal reproduction coefficient in FSTC 
are primarily determined by the Doppler effect. The thermal 
reproduction coefficient in FSTC, a//(salt), is a negative feed- 
back. As seen from Figure 5 (c), the magnitude of thermal 
reproduction coefficient in FSTC increases first and then de- 
creases as fuel salt volume fraction increases, and shows a 
monotonic increasing tendency as assembly size increases. 

The thermal utilization coefficient in FSTC is a product of 
fuel salt Doppler effect and fuel salt density effect. Because 
both the thermal utilization coefficient caused by the Doppler 
effect and that caused by the fuel salt density effect are nega- 
tive, the thermal utilization coefficient in FSTC, af (salt), is as 
well. From Figure 5 (d), at a constant assembly size, the mag- 
nitude of thermal utilization coefficient in FSTC decreases, 
owing primarily to a decrease in thermal utilization coeffi- 
cient caused by the fuel salt density effect as fuel salt volume 
fraction increases while fuel salt channel spacing decreases. 
For a constant fuel salt volume fraction, the magnitude of 
thermal utilization coefficient in FSTC increases slightly due 
to a slightly increased thermal utilization coefficient caused 
by the Doppler effect, with an increase in assembly size cor- 
responding to a larger fuel salt channel spacing. In short, the 
negative thermal utilization coefficient in FSTC is positively 
correlated with fuel salt channel spacing. The larger the fuel 
salt channel spacing, the stronger the negative thermal utiliza- 
tion coefficient in FSTC. 

Figure 5 (e) displays the change in FSTC with fuel salt 
volume fraction and assembly size. In the over-moderated re- 
gion, the feedback and the magnitude of FSTC are primarily 
determined by the thermal utilization coefficient, and FSTC 
is always negative. The magnitude of FSTC decreases as fuel 
salt volume fraction increases, and there is very little differ- 
ence in magnitudes between different assembly sizes. In the 
under-moderated region, the feedback of FSTC, mainly due 
to fast fission coefficient, resonance escape coefficient, and 
thermal utilization coefficient, could be negative or positive. 
In this case, as fuel salt volume fraction increases, the nega- 
tive FSTC decreases first due to a gradual weakening nega- 
tive thermal utilization coefficient, and then increases due to 
the enhancement rate of fast fission coefficient being faster 
than that of resonance escape coefficient. Furthermore, the 
magnitude of FSTC declines with increasing assembly size, 
owing primarily to a decreasing negative fast fission coeffi- 
cient. When the magnitude of positive resonance escape co- 
efficient is larger than the magnitude of negative fast fission 
coefficient, FSTC presents a possibility of positive feedback 
especially when assembly size is greater than or equal to 21 
cm. It can be concluded that the variation of FSTC with fuel 
salt volume fraction is primarily caused by the fuel salt den- 
sity effect, whereas the variation of FSTC with assembly size 
is primarily caused by the fuel salt Doppler effect. 


E. Moderator temperature coefficient 


As graphite moderator temperature rises, the excitation of 
target nuclei (graphite nuclei) becomes more easily activated, 
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the non-elastic scattering cross section of target nuclei be- 
comes larger, neutrons are more likely to obtain the sound 
energy of the scattered target nuclei, neutron spectrum moves 
toward the region where neutron energy is high. And a hard- 
ener neutron spectrum is revealed. The graphite moderator 
temperature coefficient (MTC) can also be divided into four 
parts, fast fission coefficient, resonance escape coefficient, 
thermal reproduction coefficient, and thermal utilization coef- 
ficient, which are caused by the graphite temperature effect. 
Figure 6 depicts the variations of MTC and its four reactiv- 
ity coefficients as a function of fuel salt volume fraction and 
assembly size. 


As graphite moderator temperature rises, neutron spectrum 
hardens, the fast neutron multiplication effect increases, and 
the fast fission factor rises. The fast fission coefficient in 
MTC, a7(gra), is positive. The magnitude of fast fission coef- 
ficient in MTC is primarily determined by the difference be- 
tween the variation in U-235’s total production and that in U- 
235’s thermal production, and is closely related to the fuel salt 
channel spacing. For a fixed assembly size, as fuel salt vol- 
ume fraction increases, fuel salt channel spacing decreases, 
neutron spectrum becomes harder and the proportion of ther- 
mal production in total production for U-235 decreases, the 
difference between the variation in U-235’s total production 
and that in U-235’s thermal production becomes larger, and 
the magnitude of fast fission coefficient in MTC increases 
(Figure 6 (a)). In the over-moderated region, for a constant 
fuel salt volume fraction, increasing assembly size results in 
a hardening neutron spectrum. In this case, the larger as- 
sembly size, the greater the parasitic thermal absorption of 
graphite, weakening the difference between the variation in 
U-235’s total production and that in U-235’s thermal produc- 
tion. And then, the magnitude of fast fission coefficient in 
MTC decreases slightly as assembly size increases. In the 
under-moderated region, when the constant fuel salt volume 
fraction increases by assembly size, fuel salt channel spacing 
increases, neutron spectrum softens, the difference between 
the variation in U-235’s total production and that in U-235’s 
thermal production diminishes, and the fast fission coefficient 
in MTC weakens. In brief, the magnitude of fast fission co- 
efficient in MTC is inversely proportional to fuel salt channel 
spacing. 

A hardening neutron spectrum caused by an increasing 
graphite temperature results in a low microscopic absorption 
cross section of graphite and a negative change in thermal 
absorption of graphite. The change in thermal absorption of 
fuel salt is affected by the fuel salt channel spacing. When 
fuel salt channel spacing is small, a hardener neutron spec- 
trum weakens the thermal absorption of fuel salt, resulting in 
a negative change in thermal absorption of fuel salt, and then 
the sum of the thermal absorptions of fuel salt and graphite 
is still negative. When fuel salt channel spacing is large, 
the average distance traveled by neutrons between two col- 
lisions increases, allowing more neutrons to pass through the 
epi-thermal region and slowly transition to thermal neutrons. 
The change in thermal absorption of fuel salt is positive. In 
this case, the decrease of thermal absorption of graphite ex- 
ceeds the increase of thermal absorption of fuel salt due to a 
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Fig. 6. Variations of different reactivity coefficients caused by the 
graphite moderator temperature effect as function of fuel salt vol- 
ume fraction and assembly size. (a) a7(gra): fast fission coeffi- 
cient in MTC; (b) af(gra): resonance escape coefficient in MTC; 
(c) af(gra): thermal reproduction coefficient in MTC; (d) at (gra): 
thermal utilization coefficient in MTC; and (e) ar(gra): MTC. 


larger graphite volume fraction. Hence, the sum of the ther- 
mal absorptions of fuel salt and graphite remains negative, 
and a negative resonance escape coefficient in MTC, af (gra), 
is revealed for all combinations of fuel salt volume fraction 
and assembly size. It can be concluded from Figure 6 (b) 
that the magnitude of resonance escape coefficient in MTC as 
a function of fuel salt volume fraction and assembly size is 
strongly related to the fuel salt channel spacing. If fuel salt 
channel spacing decreases due to a change in fuel salt volume 
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of fuel salt gradually changes from positive to negative. It 
means that the net decrease in thermal absorption of fuel salt 
grows larger and larger. Meanwhile, the increase in thermal 
absorption of graphite grows larger as well. Thus, the sum of 
changes in thermal absorptions of fuel salt and graphite en- 
hances the resonance escape coefficient in MTC. To sum up, 
the negative resonance escape coefficient in MTC is generally 
inversely related to the fuel salt channel spacing. 


For most of the combinations of fuel salt volume frac- 
tion and assembly size, an increasing graphite temperature 
causes a hardening neutron spectrum and a decrease in both 
U-235’s thermal production and U-235’s thermal absorption. 
Although U-235’s thermal production is larger than U-235’s 
thermal absorption, the former decreases faster than the lat- 
ter. The variation in U-235’s thermal production decreases 
more than that in U-235’s thermal absorption, and the ther- 
mal reproduction factor decreases, causing a negative thermal 
reproduction coefficient in MTC (a-(gra)). However, there 
are two exceptions correspond to very large fuel salt chan- 
nel spacing and very small fuel salt channel spacing. When 
fuel salt channel spacing is very small (e.g., L= 3 cm and 
VF=40%), a harder neutron spectrum causes a decrease in 
U-235’s thermal production and U-235’s thermal absorption, 
but because graphite is relatively thin, the probability of ther- 
mal neutron generation is reduced. The values of U-235’s 
thermal production and U-235’s thermal absorption are close, 
and the difference between them exceeds the difference be- 
tween the decreases of U-235’s thermal production and that 
of U-235’s thermal absorption. And then, the variation in U- 
235’s thermal production decreases less than that in U-235’s 
thermal absorption, and thus the thermal reproduction factor 
decreases and a positive thermal reproduction coefficient in 
MTC is presented. For extremely large fuel salt channel spac- 
ing (e.g., L=24 cm and VF=1%), neutrons passing through 
thick graphite between two collisions with nuclear fuel have 
a higher chance of slowing down to thermal neutrons. U- 
235’s thermal production and U-235’s thermal absorption are 
increasing, but the gap between them is wider than the gap be- 
tween the increment of U-235’s thermal production and that 
of U-235’s thermal absorption. The variation in U-235’s ther- 
mal production increases less than that in U-235’s thermal 
absorption, the thermal reproduction factor decreases, and a 
negative thermal reproduction coefficient in MTC emerges. 
According to the findings of MTC [12], the magnitude of ther- 
mal reproduction coefficient in MTC is inversely related to 
the U-235 enrichment and the EALF, mainly due to the com- 
petition between the change of thermal production and that 
of thermal absorption. On the one hand, as shown in Figure 
2 (c) and Figure 2 (d), the required critical enrichment of U- 
235 decreases first and then increases while EALF increases 
monotonously with increasing fuel salt volume fraction at a 
fixed assembly size. After superimposing the two effects, the 
magnitude of thermal reproduction coefficient in MTC with 
increasing fuel salt volume fraction for a constant assembly 
size is as follows (Figure 6 (c)): increases slowly at first, pri- 
marily due to a decreasing U-235 enrichment, and then de- 
creases rapidly due to an increasing U-235 enrichment and an 
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increasing EALF, and eventually turns to a positive feedback 
when fuel salt channel spacing is very small. The correlation 
between thermal reproduction coefficient in MTC and assem- 
bly size, on the other hand, is related to the moderated re- 
gions. As assembly gets closer to the over-moderated region, 
both U-235 enrichment and EALF increase as assembly size 
increases, and the magnitude of thermal reproduction coeffi- 
cient in MTC decreases. Also, as assembly gets closer to the 
under-moderated region, the values of U-235 enrichment and 
EALF decrease as assembly size increases, and the magnitude 
of thermal reproduction coefficient in MTC increases. 


As the temperature of graphite moderator rises, neutron 
spectrum hardens, causing thermal absorption of graphite to 
fall. The fuel salt channel spacing influences the change in 
thermal absorption of fuel salt. The neutron spectrum effect 
dominates for a smaller assembly size (e.g., L=3 cm). The 
harder neutron spectrum reduces the thermal absorption of 
fuel salt, including the thermal absorptions of heavy metal 
and carrier salt. Because the graphite thickness for neutrons 
passing through between two adjacent fuel salt channels is 
relatively thin in this case, the decrease of thermal absorption 
of heavy metal is faster than that of the sum of the thermal 
absorptions of carrier salt and graphite, and the decrement of 
the variation of thermal absorption of heavy metal is higher 
than that of the sum of the thermal absorptions of carrier 
salt and graphite, and then the thermal utilization factor de- 
creases, and thus the thermal utilization coefficient in MTC, 
a$ (gra), shows a negative feedback. For a larger assembly 
size (e.g., L=24 cm), the feedback of thermal utilization co- 
efficient in MTC is divided into two cases according to the 
moderated regions. On the one hand, when assembly tends 
to the over-moderated region (e.g., L=24 cm and VF=1%), 
the probability of a neutron slowing to a thermal neutron in- 
creases, as does the thermal absorption of fuel salt, increas- 
ing both thermal absorption of heavy metal and that of carrier 
salt. Because of the relatively thicker graphite moderator, the 
decrease of thermal absorption of graphite is greater than the 
increase of thermal absorption of carrier salt in this case. The 
sum of thermal absorptions of carrier salt and graphite has a 
decreasing trend, and then the increment of the variation of 
thermal absorption of heavy metal and the decrement of the 
variation of the sum of the thermal absorptions of carrier salt 
and graphite make thermal utilization factor increases, and 
thus the thermal utilization coefficient in MTC is positive. On 
the other hand, when assembly tends to the under-moderated 
region (e.g., L=24 cm and VF=40%), a harder neutron spec- 
trum decreases the thermal absorption of heavy metal and the 
sum of the thermal absorptions of carrier salt and graphite. 
In this case, since the thermal absorption of heavy metal is 
much larger than the sum of the thermal absorptions of car- 
rier salt and graphite, the variation of thermal absorption of 
heavy metal is larger than the variation of thermal absorp- 
tion of carrier salt and graphite, which causes a positive ther- 
mal utilization coefficient in MTC. Therefore, the thermal uti- 
lization coefficient in MTC exhibits a positive feedback for a 
larger fuel salt channel spacing or a negative feedback for a 
smaller fuel salt channel spacing. It can be seen from Figure 
6 (d), as fuel salt volume fraction increases for a fixed assem- 


oss bly size, the magnitude of thermal utilization coefficient in 
s9 MTC decreases because Maxwell spectrum is no longer visi- 
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ble. At a constant fuel salt volume fraction, as assembly size 
increases, the fuel salt channel spacing increases, the graphite 
thickness through which fast neutrons released from fuel salt 
increases, and the change in thermal absorption of graphite 
becomes more visible. The ratio of the sum of the thermal 
absorptions of carrier salt and graphite to the thermal absorp- 
tions of all materials is less than 1.0 but grows larger, and 
then the thermal utilization coefficient in MTC changes from 
less negative to more positive. When combined with the pre- 
ceding analysis, the thermal utilization coefficient in MTC is 
more negative for a smaller fuel salt channel spacing; other- 
wise, it is less negative or even positive. 

Figure 6 (e) depicts the change in MTC with fuel salt vol- 
ume fraction and assembly size. Because the fast fission co- 
efficient and the resonance escape coefficient have opposite 
feedbacks but similar values, the feedback and the magnitude 
of MTC are primarily determined by the thermal reproduc- 
tion coefficient and the thermal utilization coefficient. The 
feedback of MTC could be negative for an assembly with a 
smaller fuel salt channel spacing or positive for an assem- 
bly with a larger fuel salt channel spacing. The variation of 
MTC with fuel salt volume fraction can be divided into two 
regions: L<12 cm and L>15 cm. When assembly size is less 
than or equal to 12 cm, because the magnitude of thermal re- 
production coefficient is larger than thermal utilization coef- 
ficient, MTC is negative. With an increase in fuel salt volume 
fraction, the magnitude of MTC first increases due to an in- 
crease in thermal reproduction coefficient, and then decreases 
due to the fact that the decreasing rate of thermal reproduc- 
tion coefficient is faster than that of thermal utilization coef- 
ficient. When assembly size is greater than or equal to 15 cm, 
MTC changes from a decreasing positive feedback effect to 
an increasing negative feedback due to a decreasing positive 
thermal utilization coefficient, and then to a decreasing nega- 
tive feedback due to a decreasing thermal reproduction coeffi- 
cient. Additionally, the difference in MTC between different 
assembly sizes is primarily caused by the thermal utilization 
coefficient. When fuel salt volume fraction is less than or 
equal to 25%, the negative MTC weakens and gradually turns 
into stronger positive feedback as the positive thermal utiliza- 
tion coefficient grows stronger. When fuel salt volume frac- 
tion is larger than 25%, the negative MTC increases slightly 
with increasing assembly size, owing to a slightly enhanced 
thermal reproduction coefficient. Overall, the larger fuel salt 
channel spacing, the more likely MTC is to be positive. 


F. Temperature coefficient of reactivity 


Similarly, when the temperatures of fuel salt and graphite 
in the assembly change, fast fission coefficient, resonance es- 
cape coefficient, thermal reproduction coefficient, and ther- 
mal utilization coefficient are contributed to the total temper- 
ature coefficient of reactivity. Their variations with fuel salt 
volume fraction and assembly size are shown in Figure 7. 

The fast fission coefficient in TCR, a7(tot), shown in 


12 


== L=12 cm - L=21 cm J 
= L=15 cm —— L=24 cm 


10 


(tot) (pem/K) 


g 
T 


(or 
1 

N 

(=) 


20 30 


O Aà œ 


af(tot) (pcm/K) 


a(tot) (pem/K) 


af (tot) (pcm/K) 


a,(tot) (pem/K) 


VF (%) 


Fig. 7. Variations of different reactivity coefficients caused by the to- 
tal temperature effect as function of fuel salt volume fraction and as- 
sembly size. (a) af (tot): fast fission coefficient in TCR; (b) af (tot): 
resonance escape coefficient in TCR; (c) ar (tot): thermal reproduc- 


tion coefficient in TCR; (d) af (tot): thermal utilization coefficient 
in TCR; and (e) œr(tot): TCR. 


103 Figure 7 (a) is a sum of fuel salt temperature effect and 
104 graphite temperature effect. It can be divided into two regions 
105 by changing fuel salt volume fraction and assembly size, 
1006 corresponding to the over-moderated region and the under- 
1007 Moderated region. In the over-moderated region, the positive 
1008 fast fission coefficient in MTC is stronger than the negative 
1009 fast fission coefficient in FSTC, the fast fission coefficient in 
1010 TCR is positive. Due to a change in fast fission coefficient 
1011 in MTC, the positive fast fission coefficient in TCR decreases 
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with increasing fuel salt volume fraction or increasing assem- 
bly size. In the under-moderated region, the negative fast fis- 
sion coefficient in FSTC is stronger than the positive fast fis- 
sion coefficient in MTC, the fast fission coefficient in TCR is 
negative and strengthes with increasing fuel salt volume frac- 
tion or decreasing assembly size. 


The resonance escape coefficient in TCR, af (tot), shown in 
Figure 7 (b) is also a sum of fuel salt temperature effect and 
graphite temperature effect. A combination of a smaller as- 
sembly size and a smaller fuel salt volume fraction increases 
the likelihood of resonance escape coefficient in TCR being 
negative, owing to a negative resonance escape coefficient in 
MTC. With a fixed assembly size and increasing fuel salt 
volume fraction, the negative resonance escape coefficient 
in TCR decreases gradually and then becomes an increasing 
positive feedback. When fuel salt volume fraction is less than 
25%, the negative resonance escape coefficient in TCR grad- 
ually decreases and becomes an increasing positive feedback 
as assembly size increases. When fuel salt volume fraction is 
larger or equal to 25%, the positive resonance escape coeffi- 
cient in TCR decreases with increasing assembly size. 


The thermal reproduction coefficient in TCR, a-i(tot), 
shown in Figure 7 (c) is a sum of thermal reproduction coeffi- 
cients in FSTC and MTC. Because the thermal reproduction 
coefficient in FSTC is much weaker than that in MTC, the 
feedback and the magnitude of thermal reproduction coeffi- 
cient in TCR are primarily determined by the latter. The feed- 
back of thermal reproduction coefficient in TCR is negative. 
The variation of thermal reproduction coefficient in TCR with 
fuel salt volume fraction and assembly size is similar to that 
in MTC. When fuel salt volume fraction increases at a con- 
stant assembly size, the negative thermal reproduction coeffi- 
cient in TCR increases first and then decreases. When assem- 
bly size increases for a constant fuel salt volume fraction, the 
thermal reproduction coefficient in TCR either weakens in the 
over-moderated region or strengthes in the under-moderated 
region. 

As a sum of thermal utilization coefficients in FSTC and 
MTC, the thermal utilization coefficient in TCR, af (tot), 
shown in Figure 7 (d) could be negative or positive. When 
assembly size is less than or equal to 12 cm, the negative 
thermal utilization coefficient in FSTC has a greater magni- 
tude than the positive thermal utilization coefficient in MTC. 
In this case, the thermal utilization coefficient in TCR is neg- 
ative, and its magnitude decreases with increasing fuel salt 
volume fraction, owing to a decreasing thermal utilization co- 
efficient in FSTC. At the same time, the magnitude of thermal 
utilization coefficient in TCR decreases with increasing as- 
sembly size due to the fact that the thermal utilization coeffi- 
cient in MTC changes from less negative to strongly positive. 
When assembly size is greater than or equal to 15 cm, as fuel 
salt volume fraction increases for a fixed assembly size, the 
negative thermal utilization coefficient in TCR decreases and 
turns into an increasing positive feedback, owing to a decreas- 
ing negative thermal utilization coefficient in FSTC, and then 
the positive thermal utilization coefficient in TCR decreases 
due to a decreasing thermal utilization coefficient in MTC. As 
assembly size increases for a constant fuel salt volume frac- 


1070 
1071 
1072 
1073 
1074 
1075 
1076 
1077 
1078 
1079 
1080 
1081 
1082 
1083 
1084 
1085 
1086 
1087 
1088 
1089 
1090 
1091 
1092 
1093 
1094 
1095 
1096 
1097 
1098 
1099 
1100 
1101 
1102 
1103 
1104 
1105 
1106 
1107 
1108 
1109 
1110 
1111 
1112 


1113 


1114 


1115 
1116 
1117 
1118 
1119 
1120 
1121 
1122 
1123 


1124 


13 


tion, the negative thermal utilization coefficient in TCR weak- 
ens and then becomes a more strongly positive feedback. In 
a word, an assembly with a smaller fuel salt channel spac- 
ing presents higher probabilities of presenting more negative 
thermal utilization coefficient in TCR. 

Figure 7 (e) depicts the change in TCR with fuel salt vol- 
ume fraction and assembly size. The TCR is negative be- 
cause the sum of the negative fast fission coefficient, the neg- 
ative thermal reproduction coefficient, and the negative ther- 
mal utilization coefficient is greater than the positive reso- 
nance escape coefficient. As fuel salt volume fraction in- 
creases for a constant assembly size, the magnitude of TCR 
varies with fuel salt volume fraction in three cases with 5% 
and 25% as the cut-off points. When fuel salt volume fraction 
is less than or equal to 5%, the magnitude of TCR decreases 
rapidly, owing primarily to the rapidly decrease of the nega- 
tive thermal utilization coefficient, especially that caused by 
the fuel salt density effect. When fuel salt volume fraction 
is between 5% and 25%, the magnitude of TCR decreases 
slowly, owing to the fact that the change rate of the sum 
of resonance escape coefficient, thermal reproduction coef- 
ficient, and thermal utilization coefficient is faster than that 
of the negative fast fission coefficient. When fuel salt volume 
fraction is larger than 25%, the TCR strengths slowly, ow- 
ing to the fact that the growth rate of the negative fast fission 
coefficient is greater than that of the sum of the other three 
reactivity coefficients. For a constant fuel salt volume frac- 
tion, the change in TCR with assembly size is primarily due 
to a balance between fast fission coefficient and thermal uti- 
lization coefficient. As assembly size increases, the decrease 
rate of the negative fast fission coefficient is faster than that of 
the negative thermal utilization coefficient, so TCR weakens 
gradually. Finally, the magnitude of TCR decreases first and 
then increases with increasing fuel salt volume fraction, and 
it decreases monotonously with increasing assembly size. It 
should be noted that the contributions of fast fission coeffi- 
cient, resonance escape coefficient and thermal reproduction 
coefficient in TCR are primarily from fuel salt density ef- 
fect while the contribution of thermal utilization coefficient in 
TCR is mainly from fuel salt density effect and graphite tem- 
perature effect. Besides, the variation of TCR with fuel salt 
volume fraction mainly comes from the fuel salt density ef- 
fect, whereas the variation of TCR with assembly size mainly 
comes from the graphite temperature effect. 


IV. CONCLUSION 


Based on a graphite-moderated and low-enriched uranium 
fueled fuel salt assembly, the four-factor formula is used to 
evaluate the effects of geometric parameters including fuel 
salt volume ratio and assembly size on TCR. When assembly 
temperature rises, the fuel salt Doppler effect and the graphite 
moderator temperature effect harden neutron spectrum, while 
the fuel salt density effect softens it, affecting both the feed- 
back and the magnitude of TCR. 

A positive fast fission coefficient and a negative resonance 
escape coefficient are observed in fuel salt Doppler coeffi- 
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cient and graphite moderator temperature coefficient. A neg- 
ative fast fission coefficient and a positive resonance escape 
coefficient are revealed in fuel salt density coefficient. The 
feedback of thermal reproduction coefficient could be posi- 
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tive or negative in all reactivity coefficients. The feedback of 1165 


thermal utilization coefficient in Doppler coefficient and den- 
sity coefficient is negative while that in graphite moderator 
temperature coefficient may be negative or positive. Further- 
more, Maxwell spectrum, resonant neutron spectrum, and full 
neutron spectrum mainly affect graphite temperature effect, 
fuel salt Doppler effect, and fuel salt density effect, respec- 
tively. The full neutron spectrum has the greatest influence 
on the changes in fast fission factor and resonance escape 
factor. From low to high, the magnitudes of fast fission co- 
efficients are az(gra), az(dop), and a;(den). Similarly, the 
magnitudes of resonance escape coefficients are also af (gra), 
a (dop), and a‘(den) from low to high. Since the thermal 
neutron spectrum has the greatest influence on thermal repro- 
duction factor, the magnitudes of thermal reproduction coef- 
ficients are a-}(den), a}(dop), and œ}(gra) from low to high. 
The thermal utilization factor is a relatively complex vari- 
able, with magnitudes ranging from low to high for af, (dop), 


af (gra), and af (den). 


FSTC is a combined effect of fuel salt Doppler coefficient 
and fuel salt density coefficient. The fuel salt Doppler co- 
efficient is always negative. The fuel salt density coefficient 
can be either negative in the over-moderated region due to 
its negative thermal utilization coefficient or positive in the 
under-moderated region due to its positive resonance escape 
coefficient. The magnitude of FSTC decreases first and then 
increases as fuel salt volume fraction increases. This is pri- 
marily due to the fact that the fuel salt density coefficient 
shifts from decreasing negative feedback to increasing pos- 
itive feedback and then decreasing positive feedback. The 
magnitude of FSTC decreases monotonously as assembly 
size increases, mainly owing to a monotonously decreasing 
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Doppler coefficient. 

MITC can be negative or positive, depending on the fuel 
salt channel spacing. The smaller the fuel salt channel spac- 
ing, the more likely MTC is to be negative, owing to a neg- 
ative thermal reproduction coefficient. Because of a gradual 
weakened thermal reproduction coefficient and thermal uti- 
lization coefficient, the magnitude of MTC decreases as fuel 
salt volume fraction increases. As assembly size increases, 
MTC transitions from weakened negative feedback to en- 
hanced positive feedback as a result of the change in thermal 
utilization coefficient, and when fuel salt volume fraction is 
relatively lower, the influence of fuel salt channel spacing on 
MTC is more visible. 

At an assembly level, the total temperature coefficient of 
reactivity is negative, regardless of whether the region is over- 
moderated or under-moderated, which meets the safety re- 
quirements of reactor operation. The magnitude of TCR in 
the under-moderated region is smaller than that in the over- 
moderated region, and the power change rate is compara- 
tively slower, which can be conducive to reactor safety. MTC 
is more likely to be negative in the under-moderated region 
especially for an assembly with a smaller fuel salt channel 
spacing. As a result, an assembly with a smaller fuel salt 
channel spacing in the under-moderated region is preferred 
to achieve a reasonable negative TCR as well as a negative 
MTC. The thermal-hydraulic analysis and the transient anal- 
ysis are also required to assess the reasonableness of the mag- 
nitude of TCR to ensure the intrinsic safety of a liquid-fueled 
MSR. Future work will include, but is not limited to, the fol- 
lowing items: friction pressure drop, heat conduction ability 
of graphite, consequences of abnormal reactivity introduced 
into, and blockage of fuel salt. In addition to geometrical 
parameters, fuel salt composition is another crucial factor in- 
fluencing TCR for a thermal MSR, especially when fuel salt 
reprocessing and online refueling are equipped. To fully un- 
derstand the mechanism of TCR for thermal MSR, more re- 
search is needed to address how TCR is affected by both ge- 
ometric parameters and fuel salt composition. 
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